On the quark propagator singularity 



V. Sauli 1 

1 CFTP and Departamento de Fisica, Instituto Superior Tecnico, Av. Rovisco Pais, 1049-001 Lisbon, Portugal 

Using the method of Fukuda and Kugo [l|] the continuation of Euclidean solution is performed 
to the timelike axis of fourmomenta. It is shown that assumed presence of the real simple pole in 
quark propagator is not in agreement with the solution. The simple pole disappears because of the 
discontinuity in the resulting quark mass function. 
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Theoretical analysis of the singularities and the be- 
haviour of QCD Greens functions in the timelike momen- 
tum regime is still not well understood. Natural nonper- 
turbative framework for the study of infrared properties 
of QCD Green's functions is the formalism of Schwinger- 
Dyson equations (SDEs) . For a review of recent progress 
in the whole QCD see |2| . Few possibilities of the quark 
propagator behaviour have been discussed in the litera- 
ture. As quarks are confined objects in nature one of the 
most natural expectation is the absence of a real pole 
in quark propagator. Some of the models suggest that 
the real pole can be split into the complex ones. The 
complex conjugated poles have been phenomenologically 
appreciated in various studies. Furthermore, it has been 
recently realized in [3j that the structure of light quark 
propagator is not reliably described by the Lehmann rep- 
resentation. The numerical analyzes exhibit inefficiency 
of spectral representation for correct description of QCD 
chiral symmetry breaking. It is also known that when the 
strong interaction happen to lead to the highly infrared 
singular interaction then the real quark pole naturally 
disappear 

However, the absence of pole in the quark propagator 
is not the only possible mechanism of quarks confine- 
ment and the existence/absence of the real pole is still 
undetrmined [5[. The primary objective of this paper is 
to check the presence / absence of a real quark propagator 
pole by using a simple generalization of the Fukuda-Kugo 
timelike continuation of ladder fermion SDE to the case 
of QCD, e.g. the asymptotic freedom is correctly taken 
into account through the known behaviour of running 
coupling. Using usual assumptions we describe contra- 
dicting results we obtained: The quark propagator does 
not exhibit the simple pole behaviour, instead, the mass 
function is largely discontinuous at observed real branch 
point. 

The following conventions are used: the positive vari- 
ables x, y represent the square of momenta such that 
x = p 2 for timelike momenta when p 2 > 0, while 
x = —p 2 for p 2 in spacelike region. Our metric is 
9nv = diag(l, — 1, — 1, — 1). For purpose of clarity the 
mass function M is labeled as M s in the spacelike region 
of fourmomenta and as M t when evaluated for timelike 
fourmomenta (i.e., M(p 2 ) = M s {x)0{-p 2 ) + M t {x)6{jp 2 )). 

The quark propagator S can be conventionally charac- 
terized by two independent scalars, the mass function M 



and renormalization wave function Z such that 
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noting the bare fermion propagator is So = 
where SDE for the inverse of S reads 
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where simply M = B/A, A = l/Z and T is the full 
quark gluon vertex, G a @ represents gluon propagator, 
both Greens functions satisfy their own SDE. To make 
the later continuation to the timelike regime more easily 
tractable we intruduce approximations that we believe 
do not change qualitatively feature of solution. We will 
work in Landau gauge and take A = 1. The impor- 
tance of (un)presence of A can be estimated from similar 
Euclidean studies. The so called Analytic Running Cou- 
pling (ARC) 0, 0] is used to properly include the running 
of the QCD coupling. Noting that the effect of dynam- 
ical chiral symmetry breaking has been already studied 
Q within ARC-SDE combined framework. 

In dressed one loop approximation the following pre- 
scription of the SDE kernel is used 
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with ARC written via dispersion relation as 
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The correct one loop QCD running coupling at asymp- 
totically large —q 2 is ensured when p g ~ p\ 



Pi (A, Aqcd) = 
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while the infrared behaviour of ARC is modeled through 
the following modification of the spectral function: 



Pg (X) = (l + f IR (X)) Pl (X) 



fm(X) 



^ -^-j? — - 
Ce a «cd 



(0.6) 



2 



where pi(A) is given by (|0.5|) . Nonperturbative extra 
term fm(X) is chosen such that it does not affect ultravi- 
olet asymptotics. Simultaneously it leads to the expected 
QCD scaling of infrared up and down quark masses 
M Ut d(0) — Aqcd, where C — 15 numerically. In order to 
define the kernel entirely we fix Aqcd — 200MeV and 
47r//3 = 1.396, noting the later value corresponds with 
one loop ARC a(0) calculated for three active quarks. 

After making the trace, SDE (|0.2|) for the mass func- 
tion reads 



B(p 2 ) = m + 3.CU J 
B(q 2 
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where Ca = T a T a = 4/3 for SU(3) group generators. 

Using the integral representation (|0 .4|) and the follow- 
ing simple algebra 
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1 



A 



(0.8) 



we get after the Wick rotation and angular integration 
the resulting equation for B: 



B s (x) = to + 
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solution solely. Within presented model the mass func- 
tion of the light and and heavy quarks has been calcu- 
lated. The current bare masses were chosen to be m u> d = 
5 MeV, m s = 100 Me V, m c = \GeV 7 m b = 3.5 GeV. 
The resulting mass functions are plotted in Fig. [TJ 

To make a continuation of Euclidean Greens function 
to the physical timelike momenta is a rather cumber- 
some task. The assumptions are indispensable and not 
rarely if checked consequently, they appear not justified 
in many physically important cases. In what follows we 
will show that the assumption of Fukuda-Kugo continu- 
ation in QCD is another case. 

Changing trivially the order of the integrations we can 
immediately follow the Fukuda-Kugo prescription and 




The gap Eq. (|0.9[) contains potential UV divergence 
and can be renormalized after a suitable regularization. 
Since low relevance for our topic, the explanation is sim- 
plified and the renormalization issues are not discussed. 
In presented work we do not renormalize at all and solve 
the regularized gap equation. For this purpose we fol- 
low hard cutoff scheme and introduce the upper integral 
regulator Ah, the numerical data has been obtained for 
A# = c 8 Aqcd- Let us mention that we have explicitly 
checked that renormalization has only marginal numer- 
ical effect when comparing to the properly regularized 



FIG. 1: Mass functions of the quarks. The thin lines are the 
results obtained by the integration through the absorptive 
part of M in assumed dispersion relation as described in the 
text. 



write down the result for the timelike momentum. Con- 
tinued SDE for the positive timelike square of fourmo- 
menta can be written in the following way: 



B t (x) = B{x)-I{x)+i 
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FIG. 2: U,D,S quark mass functions for positive p 2 . Thick 
(thin) lines represent the real (imaginary) parts. The dotted 
line is liner function in p. 

where the function K(a,b, c) is defined in Eq. (|0.10p . 
to is a pole mass, and the function X is defined as 

X( X ;y, Z ) = ^lE^EM. (0.12) 

x 

Eq. (|0.10p represents the integral equation for the mass 
function B t {x) defined at the timelike axis x. It consists 
of the dominant term B(x) represented by regular inte- 
gral which including the mass spacelike regime defined 
function B s . It also contains the principal value integral 
that includes the complex function B t itself. As usu- 
ally, absorptive part of B t is generated when crossing the 
branch point. 

The gap equation has been solved by the standard nu- 
merical iteration. To check a numerical stability several 
integrators have been used to perform the integration nu- 
merically, e.g. the Gaussian and Simpson ones. The ob- 
served form of the discontinuity in numerically obtained 
quark mass function strongly indicates that the singular- 
ity of quark propagator is softened such that 

Um p 2^ m i(p 2 — m 2 )S(p) — 0, (0.13) 

which is in contradiction with the assumption. The limit 
(10.13)) is well defined since the left and right limits coin- 
cide. Note taht this is not a case of S because of discon- 
tinuity in M t . The function M t becomes complex above 
the to with ImM starting from zero at this point. How- 
ever, the right limit, if considered for the real part of M, 
is a bit model dependent, however we should stress that 
it has never leaded to the singular propagator. There 
are two possible realizations observed numerically which 
spectacularly differ by the behaviour of M t in the right 
vicinity of to. Using a large UV cutoff one can find that 



FIG. 3: C,B and S quark mass functions for positive p 2 . 
Thick (thin) lines represent the real (imaginary) parts. The 
dotted line is liner function in p. 

the function Mt cuts the linear function p for some p > to 
but ImMt is always nonzero at that point. The second 
possibility that the inequaility p — M t (p) ^ is hold for 
any p is observed for low scale cutoff. Such situation 
is displaed in Fig. 4 for the case A = 5GeV and light 
quark mass to = 5MeV. Of course, the later case is not 
real QCD, but this is rather qualitative test of cutoff in- 
dependence. As mentioned the left limit of M t is less 
affected by the model changes, here it has been lowered 
about 50MeV. Thus in any case, one can conclude that 
the shape and perhaps complexity of the function M t in 
the vicinity of its discontinuity prohibit the appearance 
of the real pole in the dressed quark propagator. 

Termed in other way, we obtain finite quark propagator 
with discontinuous mass function at the point where the 
ordinary threshold would be expected. If one need, the 
discontinuity can be classified 

Um p 2^ m 2_TrS(p) ~ (0.14) 

where A = B t (m^) — to is the "left discontinuity" of 
the function B t at the point to. The point to is just 
defined by this discontinuity. Its numerical value has 
been searched by minimization of A during the itera- 
tion procedure for each quark flavor separately. The nu- 
merical values we have found start from approximately 
A u <i ~ lOOMeV for vanishing current masses and slowly 
grows up to be A Cj {, ~ 400AfeU for heavy quarks. The 
ratio A/to decreases for heavier flavors, which is agree- 
ment with expected supression of nonperturbative effects 
in the case of heavy quarks. The numerical solutions are 
displayed in Fig. [5] for the light quarks u,d,s and in 
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Fig. [3] for the heavy flavors, in the later case, the mass 
function for the strange quark is added for comparison. 

In a strict sense the observed solutions should not be 
taken seriously, since based on wrong assumptions. On 
the other hand, the timelike continuation up to the dis- 
continuity point m should be trustworthed. Then we can 
conludc that there is strong indication of disappearance 
of pure real pole at all, no matter what the continuation 
above the point m is. The observation of mass gap A is 
in agreement with the confinement of quarks. Hadrons 
can never dissociate into the free quarks. A development 
of a new techniques for more complete study of timclike 
behaviour of QCD Greens function are beeing looked for. 
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FIG. 4: The detailed picture of discontinuity, the details are 
described in the text. 
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